Bidiagonalization with Parallel Tiled Algorithms

Block and Parallel Versions of One-Sided Bidiagonalization

Two new algorithms for one-sided bidiagonalization are presented. The first is a block version which improves execution time by improving cache utilization from the use of BLAS 2.5 operations and more BLAS 3 operations. The second is adapted to parallel computation. When incorporated into singular value decomposition software, the second algorithm is faster than the corresponding ScaLAPACK rout...

Adaptive Parallel Householder Bidiagonalization

With the increasing use of high-resolution multimedia streams and large image and video archives in many of today’s research and application areas, there is a growing need for multimedia-oriented highperformance computing. As a consequence, a need for algorithms, methodologies, and tools that can serve as support in the (automatic) parallelization of multimedia applications is rapidly emerging....

Parallel Computation of Spectral Portraits of Matrices by Bidiagonalization

Parallel Band Two-Sided Matrix Bidiagonalization for Multicore Architectures

The objective of this paper is to extend, in the context of multicore architectures, the concepts of algorithms-by-tiles [Buttari et al., 2007] for Cholesky, LU, QR factorizations to the family of twosided factorizations. In particular, the bidiagonal reduction of a general, dense matrix is very often used as a pre-processing step for calculating the singular value decomposition. Furthermore, i...

Implementation of Parallel Numerical Algorithms Using Hierarchically Tiled Arrays

In this paper, we describe our experience in writing parallel numerical algorithms using Hierarchically Tiled Arrays (HTAs). HTAs are classes of objects that encapsulate parallelism. HTAs allow the construction of single-threaded parallel programs where a master process distributes tasks to be executed by a collection of servers holding the components (tiles) of the HTAs. The tiled and recursiv...